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ABSTRACT 

I develop a method for assessing the ability of an instrument, coupled with 
an observing strategy, to measure the angular power spectrum of the cosmic 
microwave background (CMB). It allows for efficient calculation of expected 
parameter uncertainties. Related to this method is a means of graphically 
presenting, via the "eigenmode window function", the sensitivity of an ob- 
servation to different regions of the spectrum, which is a generalization of 
the traditional practice of presenting the trace of the window function. I 
apply these techniques to a balloon-borne bolometric instrument to be flown 
this spring (MSAM2). I find that a smoothly scanning secondary is better 
than a chopping one and that, in this case, a very simple analytic formula 
provides a good (40% or better) approximation to expected power spectrum 
uncertainties. 



1 Introduction 



Precision measurement of the angular power spectrum of the cosmic mi- 
crowave background (CMB) promises enormous scientific returns. The ex- 
pected results for the next generation of satellite experiments have been 
studied extensively [p], [|. If one assumes that structure was formed by 
gravitational instability then MAP g] and COBRAS/SAMBA @ should 
both measure Q, the Hubble constant and other cosmological parameters to 
better than a few per cent. 

These large sky coverage, map-making observations lend themselves to 
easy analytic evaluation. However, there are very important balloon-borne 
and ground-based observations to be done over the next few years for which 
the determination of expected power spectrum uncertainties is not as straight- 
forward. The method presented here is in principal analytic but the necessary 
high-dimensional linear algebra requires numerical work in practice. 

One influence on the ability of any instrument to measure the microwave 
background is its ability to separate the CMB from confusing astrophysical 
foregrounds. Much work has been done along these lines || [|. In contrast, 
relatively little attention has been paid to the choice of spatial observing 
strategy and its effect on expected power spectrum uncertainty. The benefits 
of high angular resolution have been emphasized |]J and it is well-known that 
a rough guide to optimal sky coverage in a fixed time is that which gives a 
signal-to-noise ratio per independent pixel of unity. However, the pixels are 
not independent, necessitating a more sophisticated analysis. 

One reason for this discrepancy between attention paid to frequency 
strategy versus spatial strategy is that while the spectrum of the CMB and 
confusing foregrounds have been known for a long time (at least roughly), it is 
only recently we have gained some knowledge of the angular power spectrum 
across a large range of angular scales. Prior to any detection of anisotropy, 
the only possible guidance on spatial strategy was that the optimal number 
of points to observe is thirteen 0. Now that our knowledge has improved 
substantially, it should be used for the choice of observing strategies. 

In section two I present the method and discuss the virtues of the eigen- 
mode window function. In section three I derive a simple analytic approxi- 
mation for power spectrum uncertainty. In section four I apply the method 
to the case of MSAM2 as an example. I use the example to make a quantita- 
tive argument for a smoothly scanning secondary as opposed to a chopping 
secondary, to present the eigenmode window function, and to compare the 
calculated power spectrum and cosmological parameter uncertainties with 
those from a simple analytic formula. An appendix details the calculation of 
the theory and noise covariance matrices and also suggests a useful prescrip- 
tion for normalizing window functions. 
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2 Semi-analytic Method 

The data set can be modeled as consisting of signal and noise, Di = Si + iVj. 
Its covariance is: 



(DtDj) = (SiSj) + (NiNj) (1) 
-is , r<N 



+ C» (2) 



For the simple case where the are measurements of the temperature in 
directions enumerated by i, the signal covariance matrix is only a function 
of the angle separating i and j, 9 if 

= C(%) = E 2 -^C l P l {cos9 ij ) (3) 

where (7/ is called the angular power spectrum. For reasons to be discussed 
below, the data are generally not direct measurements of the temperature 
in a given direction, but are linear combinations of temperatures in several 
directions. Calculation of C s and C N in the general case is discussed in the 
appendix. 

In general, the signal covariance matrix depends only on the observ- 
ing strategy and the angular power spectrum G\. We can think of the 
angular-power spectrum as a function of parameters, a p , either physical 
(flo, us-, h, etc.) or phenomenological. The task at hand is to estimate, for 
a given experiment, the uncertainties to expect on the parameters a p . 

Our task is simplified by transforming the data so that both the signal 
and noise covariance matrices are diagonal. The desired transformation is 
the Karhunen-Loeve [|| transformation or, as it has become known in CMB 
phenomenology, the signal-to- noise eigenmode transformation JPJ. Bond \TU 
and White and Bunn [|TT|] used it in their analyses of the COBEDMR maps. 
Vogeley and Szalay have shown how it is useful for analysis of galaxy redshift 
surveys |T2j . The transformation is a non- unitary mapping of the data in the 
pixel basis into the "signal-to-noise eigenmode" basis. The transformation 
to D[ (the data in the S/N basis) is as follows: 

D[ = R ti ((C N )- 1 '*) jk D k (4) 

where C^j = (NiNj), R is the rotation matrix that diagonalizes 

M = (C N )- 1 / 2 C S (C N )- 1 ^ 2 and = (SiSj). It is straightforward to show 

that 

(D^D') = (Aj + 1)^ (5) 
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where Aj are the eigenvalues of the S/N matrix, M. Thus we have the 
desired transformation, since the signal and noise covariance matrices are 
now diagonal. 

We can now think of the experiments as measuring the Aj which are a 
function of the power-spectrum and therefore of the parameters, a p . We can 
build a quadratic estimator for Aj from D[ by squaring it and subtracting 
unity: 

Af = Df - 1 (6) 

This estimator has variance 2(Aj + l) 2 since D[ is a Gaussian random variable 
with variance Aj + 1. 

From here it is straightforward to calculate uncertainties on the parame- 
ters. First calculate the curvature matrix a: 



dX, OX 



a: 



y (7) 



^ cr 2 (X, ) da p da 



where cr 2 (Aj) = 2(Aj+l) 2 . Then invert to get the desired result, the parameter 
covariance matrix [13|, |T 



C^a^L. (8) 



pp \ / PP > 

Equation ||| is strictly true only if our estimates of dp 3X6 Gaussian random 
variables which will be the case if the eigenvalues depend linearly on a p . In 
general the Aj have non-linear dependences on physical parameters, but if 
the data constrain the parameters sufficiently then the non-linear effects will 
be unimportant. 

The covariance matrix is simply given by the inverse of the curvature 
matrix for the case when we have no prior information on the parameters. 
At the other extreme, if the parameters other than a p are perfectly known 
then a 2 (a p ) = l/a pp . The intermediate case is easily treated if the prior 
information is expressed as a covariance matrix Cp r i or - Then 



a 2 1 



a p ) = ((a + Cj\ ') . (9) 

To summarize, the procedure is straightforward. First choose a 
parametrized theory and calculate the signal and noise covariance matrices. 
This step requires a specific choice of theory parameters. Calculate the rota- 
tion matrix, R, the eigenvalue spectrum, Aj, and its derivatives. Calculate the 
curvature matrix, add any prior information and invert to get the estimated 
parameter covariance matrix. 

The imaginary analysis of the data that I used above to derive Eq. |7| 
neglected off-diagonal terms that would be included in a full likelihood 



3 



analysis. A Taylor series expansion in the log of the likelihood function, 
about its maximum, leads to the result [|1^ 

2a pp , = tr[(C s + C»)-^{C S + C N )-^\ (10) 
where tr stands for the trace. With a little algebra this becomes: 



da, n l 



a 



VP 



h3 



2(A i + l)(A i + l 



The difference between the two results (Eq. [8] and Eq. O ) is simple to 
understand. The imaginary analysis above neglected the offdiagonal terms 
that arise when we move in parameter space, off of the parameters for which 
the signal matrix is diagonalized. In other words, although A is diagonal by 
design, its derivatives are not. Since the maximum likelihood is the "best" 



[ |16|j estimator, the parameter variances following from Eq. Ill] will be smaller 
than those from Eq. & I will refer to the use of Eq. 18 as the diagonal 



approximation and to Eq. 11 as "exact", with quotations due to the implicit 



approximation of linear parameter dependence already mentioned. 



Eq. TT| is an improvement on Eq. |10| because of the compression it allows; 
the modes below some level of signal-to-noise can be discarded. Because it 
neglects the off-diagonal terms, the diagonal approximation is even more 
convenient and in the applications that follow I show that the information 
loss from neglecting the off-diagonal terms is small. 

Another advantage of the accuracy of the diagonal approximation is that 
it allows for easy visualization of how the observations will probe the spec- 
trum. The window function Wuj relates the power spectrum and its deriva- 
tives to the signal covariance matrix and its derivatives: 

da p ^ da p Air ^ 

In the pixel basis the window function is complicated. The usual practice 
for indicating what region of the spectrum is probed is to plot only its trace 
- a procedure which neglects possibly important off-diagonal terms. It is 
possible to do much better in the eigenmode basis. From Eq. [E2] follows: 

da p -^da p H > 4vr {U} 

where the prime indicates the eigenmode basis. In this basis one can safely 
ignore the off-diagonal elements, reducing the number of relevant dimensions 
to two (one mulitpole moment index plus one pixel index), allowing for 
visualization of all the important parts of the window function. 
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3 Simple Analytic Methods 



For a map with uniform full-sky coverage and Gaussian angular resolution, 
<Jb, the eigenmodes are the spherical harmonics with eigenvalues equal to 
wC\e~ a b where w is the weight per solid angle. There are 21 + 1 modes for 
each I and thus cr 2 (Aj) = 2(A« + l) 2 leads to the formula (derived by other 
means in fllfl): 



^^ — (Q + w-'e 1 ^) . (14) 

The effect of observing only a fraction of the sky, /, can be approximately 
described by increasing the variance of C\ by f -1 since the number of modes 
is roughly proportional to the area [§, O, |18[ - Jungman et al. @| used this 



formula and the analogs of Eq. [5] and |8pto calculate standard errors for an 
eleven parameter adiabatic gravitational instability model. The partial-sky 
corrected version of Eq. [14| must be used with care, however. The modes 
are no longer spherical harmonics and therefore any estimate for Ci will be 
correlated with that of C[ - a correlation with range Al m 2ir/6 where 9 is 
a linear dimension of the observed field ||19|| . The equation only makes sense 
when the spectrum is binned with Al > 2%/ 9. 



4 Application to MSAM2 

In this section I apply the above formalism to the particular case of the second 
Medium Scale Anisotropy Measurement instrument (MSAM2). The MSAM 
instruments are balloon-borne off-axis Cassegraine telescopes with bolometric 
radiometers. The MSAM1 instrument is described in pD| |. Detections have 



been reported from two of the flights [2~I| with the second flight confirming 
the results of the first [p2|. The MS AM2 instrument uses the same gondola 
as MSAM1 but has a different radiometer and secondary mirror. 

The largest component of atmospheric contamination depends only on 
elevation and is slowly varying in time. Thus a standard technique to reduce 
atmospheric contamination is to rapidly sample a stretch of sky at constant 
elevation by motion of a secondary mirror from — 9 C to 9 C . Only linear combi- 
nations of the data that have no sensitivity to a spatially homogeneous signal 
are kept for further processing. In some cases those combinations sensitive 
to a gradient are discarded as well. Each linear combination corresponds to a 
"synthesized antenna pattern" . The MSAM1 secondary motion was a three 
point chop. From the time stream, two antenna patterns were synthesized. 

The MSAM2 secondary motion will be a triangle wave pattern of period 
T. Here we model the time stream of data as follows: 
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d(t) = + £ DftOi) cos(27r/i-) + £ DJ(ft) sin(27rz/-). (15) 

where 6i is a slowly changing function of time and refers to the point on the 
sky observed when the secondary is in its central position. 

We cut the Fourier decomposition off at v — 11 — 8 because the peak-to- 
peak chop amplitude is 8 beam full-widths; higher frequency modes would 
have very little signal. Since the secondary motion is symmetric about t = 
T/2, the asymmetric Fourier modes will contribute zero signal. Thus the odd 
v components are ignored. 

If we assume that the noise in the time stream d(t) is white, then the 
noise in each of the above modes will be independent; the noise covariance 
matrix C N will be diagonal (see appendix). A better model would also have 
terms that vary in time but not in space, as is the case for instrumental and 
atmospheric drifts. Having to fit for the coefficients of these terms would 
induce correlations in the noise covariance matrix C N . Given a model of the 
drifts it is possible to calculate C N but here I ignore these effects and take the 
matrix to be diagonal. From a model of the bolometer and the foregrounds 
[§] (a one component dust model) we expect the sensitivity to CMB to be 
NETf| ~ 266 /xKv/iec. 

Although the noise matrix is diagonal, the signal matrix is not. The off- 
diagonal correlations exist for three conceptually distinct reasons. First, the 
Fourier decomposition in Eq. 15] is for functions of period T but because 
of the triangle wave motion ofthe secondary, the same stretch of sky is 
scanned twice in that one period. The D^(8i) and D s ^(8i) are the cosine and 
sine coefficients of the sky sampled from Q { — 9 C to Qi + 9 C , with fundamental 
frequency half what it would be for a Fourier decomposition. Second, because 
the central position of the chopper changes over time, the decompositions do 
not all share the same origin. Thus, even if the fundamental mode had the 
right frequency, the modes with 9i ^ 6j are not orthogonal. Third, since the 
sky is not "white" like the noise (there are intrinsic correlations) the different 
modes are correlated - even when \6i — 6j\ » 29 c . The calculation of the 
signal covariance matrix, C s , is described in the appendix. 



4.1 Optimal Motion of the Secondary 

Is it better to move the secondary in a step motion between two or three 
spots, or to smoothly scan it back and forth? Here I have addressed that 

lr The standard error of the observed temperature is given by NET/^/i b s where i b s is 
the observing time. 
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question by computing the S/N eigenvalue spectra, shown in Figure [I], for a 
three-point chop and a triangle wave. 




20 40 60 80 

i 

Figure 1: The S/N eigenvalue spectrum for a three-position secondary and 
for a smoothly scanning secondary. 

The two curves shown in Fig. 1 are both S/N spectra for observing 
strategies that are the same in all respects except for the motion of their 
secondaries. They both point at a declination of 5 = 78 degrees on the 
transit meridian as the secondary moves back and forth ±80' at constant 
declination (an approximation to the actual motion in cross-elevation). The 

beam is taken to be a Gaussian with = 20'/ (v8Zn2). The sky is observed 
in this manner for five hours. The rotation of the sky leads to coverage of a 
strip 15.6 degrees long. For a theory I took a flat spectrum with: 

51? - = (50 fiK) 2 . (16) 

The dashed curve is for the secondary that executes a three-point chop. 
The single difference and double difference signals are analyzed. The solid 
curve is for the triangle wave motion. In this case the signals from the 
12 different synthesized antenna patterns were analyzed. One can see that 
the chopping secondary is inferior to the smoothly scanning ones since its 
eigenvalue spectrum is lower for every mode. 

As a rough guide to the power-spectrum sensitivity of the experiment, 
we can simply count the number of modes with A, > 1. This is because 
modes with \ » 1 all have the same fractional error (Vz) while those with 
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Aj << 1 have very large uncertainties. For a more precise comparison we can 
perform the sum 



where A is the amplitude of the spectrum whose shape has been assumed. 
For the triangle wave spectrum <j(A) = 0.22. 

It is easy to see how sensitivity to the amplitude of the spectrum will 
change with varying sky coverage. Increasing the sky coverage by a factor 
of n will increase the number of modes by a factor of n. If the observing 
time remains fixed, then the noise in each pixel will increase, reducing each 
eigenvalue by a factor of n. From Fig. 1, we see that it would be highly 
beneficial to greatly increase the sky coverage. 

The uncertainty in the amplitude has a very shallow minimum at n = 51 
of cr{A) = 0.087. With such large sky coverage, the noise in each beam-size 
pixel is 97 /iK. This is much larger than the signal in any of the synthe- 
sized antenna patterns but is just slightly smaller than the "stare mode" 
(undifferenced) signal of 120 fiK. Thus we see that if one is solely interested 
in measuring an amplitude, setting the undifferenced signal-to-noise ratio to 
unity gives a near optimal sky coverage. However, the need to understand 
the inevitable non-idealities of the data argues against spreading the weight 
this thinly. 



4.2 Band Estimates 

The eigenvalue spectrum shows the sensitivity of the experiment to the over- 
all amplitude of the spectrum but gives no indication of sensitivity to shape. 
To understand what regions of the spectrum are being probed it is useful 
to look at the eigenmode window function, shown in Fig. |2| for the triangle 
wave observations of the previous subsection. As expected, the peaks move 
to higher values of I as the mode number increases. 

Traditionally the region probed by an experiment has been indicated by 
the diagonal elements of the pixel basis window function, Wm |f23 |. This 



is adequate if all that is important is the rms of the data. However, off- 
diagonal correlations are important as well and can significantly alter what 



regions of the spectrum are probed [24|. In particular it should be noted 



that as signal-to-noise ratio increases more information starts coming from 
the smaller scales - an effect that should be clear from Fig. [| To quantify 
how well the experiment is probing different parts of the spectrum it is useful 
to estimate the uncertainty in several bands. To simplify interpretation it is 
useful to choose the bands with widths greater than the eigenmode window 
function widths since that will reduce the correlations between estimates. 
For Fig. H we have parametrized the spectrum as 
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Figure 2: A few modes of the eigenmode window function 



(D t if 20 < I < 210 
8Ti = I D 2 if 210 < I < 400 (18) 
{D 3 if 400 < / < 810 

and taken Di = D 2 = D 3 = 50 fiK. Note that estimating the power in a given 
band is very different from estimating the "band-power" for a synthesized 
antenna pattern. In the former we are using all the data to constrain part 
of the spectrum. In the latter we are using part of the data to constrain the 
entire spectrum. 

The open triangles in Fig. |3] show the expected uncertainties for the trian- 
gle wave observing strategy described in the previous subsection. I will refer 
to this strategy as "parallel" since the sky coverage is swept out parallel to 
the motion of the secondary. The open squares are for a strategy where the 
sky is swept out perpendicular to the secondary motion. Each case has the 
same area: 5.2 sq. degrees, the same sensitivity, observing time, beam size, 
and even the same synthesized antenna patterns. The one long dimension of 
the "parallel" strategy is the reason for its superiority in the lowest spatial 
frequency bin. 

The band estimates allow us to see the effect of varying sky coverage in 
more detail. The closed symbols in Fig. |3] are the expected uncertainties 
for observing 16 patches of the sky in the same way, with the same total 
observing time. The 16 patches are assumed to be far enough apart that the 
correlations between them can be ignored. In this case, the calculation is a 
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0.1 
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Figure 3: The expected fractional standard error for the three bands for 
the "parallel" (triangles) and "perpendicular" (squares) observing strategies 
(see text). The pentagons are the corresponding results from the analytic 
approximation. Filled symbols are for sixteen times as much for the 

open symbols; observing time is fixed at 5 hours. 

simple extension of that done for the open symbols: the number of modes is 
increased by 16 and the S/N eigenvalues are all reduced by a factor of 16. 
Note that as expected, the lower the spatial frequency of the band, the more 
it benefits from extra sky coverage. 

The "parallel" and "perpendicular" strategies are two extremely different 
ways of uniformly covering 5.2 sq. degrees. Despite this fact, the analytic 
approximation (pentagons), which does not depend on how the sky is covered, 
is good for both of them to better than 40%. Where the strategies are 
most different (-Di), it splits the difference! The results from the diagonal 
approximation (not shown) differ by at most 10%. 

4.3 Cosmological Parameters 

The sky coverage planned for the MSAM2 flight is too small to allow for 
simultaneous determination of a large number of cosmological parameters. 
However, it is interesting to vary a limited set of parameters and contrast 
the results from the "exact" formula with the diagonal approximation and 
the analytic approximation. To calculate the standard errors in the Table I 
varied the quadrupole, (Q), the scalar power-law index ns and Q to t = flo+fl\ 
around their standard CDM values of (Q) = 20 fiK, ns = 1 and Q t ot = 1 
[BHj. For prior information, I assumed that COBE-DMR determined the 
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quadrupole to within 20% ||25|| . Calculations were done for nominal and 16 
times sky coverage. 



Table 1: Constraints on Cosmological Parameters 



Sky Coverage 


Method 


A(Q) 


An s 


Afi tot 


1 


Diagonal 


4.0 


0.27 


0.35 


1 


"Exact" 


4.0 


0.22 


0.28 


1 


Analytic 


4.0 


0.18 


0.22 


16 


Diagonal 


4.0 


0.15 


0.17 


16 


"Exact" 


3.9 


0.12 


0.14 


16 


Analytic 


3.9 


0.11 


0.12 



It is not surprising that the diagonal approximation does worse here than 
for the band estimates since the diagonal is equal to the "exact" when the 
parameter is a normalizing constant. But even with these cosmological 
parameters it is still good to 25%, bolstering the claim that the diagonal 
elements of the eigenmode window function hold most of the information. 
The analytic approximation, off by 20% at the worst, is even better here 
than for the band estimates. Note that the constraint on (Q) comes almost 
entirely from the prior information. 



5 Discussion 

The simple analytic approximation worked quite well for the cases stud- 
ied. The approximation will be worse for observations where the differencing 
throws away more information, the off-diagonal noise is more important, or 
the sky coverage is less uniform. Differencing will throw away more infor- 
mation when the secondary throw to beam fwhm ratio is smaller or, as was 
explicitly shown, if the secondary chops discontinuously. On the other hand, 
we can expect the analytic formula to be a very good approximation for 
partial sky map-making observations, such as those to come from Long Du- 
ration Ballooning. 

Besides a method for estimating parameter uncertainties I have also 
presented a useful tool for understanding how an observation probes the 
spectrum. Combined with a graph of the eigenvalue spectrum, the diagonal 
elements of the eigenmode window function allow one to immediately see the 
£-space resolution, the sensitivy to different regions of the spectrum, and the 
dependence of that sensitivity to varying sky coverage. The traditional use of 
the diagonal elements of the pixel basis window function can be misleading 
if the experiment has a high signal-to-noise ratio. 



11 



The bands in Fig. |3] were chosen with large widths to reduce correla- 
tions. However, these large widths mean some resolution has been discarded 
and thus this is not a good solution to the problem of presenting data in a 
compact, yet complete, manner. Further work along these lines is needed. 
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A Covariance Matrix Calculation 

Given an observing strategy and a theory, we can construct the theoretical 
correlation matrix. Here I show how that calculation is done in general. For 
simplicity of notation, we assume that the telescope is pointed at 9i for only 
one cycle of the secondary before moving on to 6j. In reality, many cycles are 
completed before the telescope's pointing has changed significantly. Assume 
that the detectors are sampled n times in one period T and di a is the direction 
on the sky observed on the a th sample of the i th secondary cycle. The multiple 
antenna patterns are synthesized by weighting the n samples of each cycle 
by the weight vector w^ a : 

1 n 

n o=l 

For example, the weight vectors for D c {6i) in Eq. |l| are 

w^a = cos {2tt ixa/n) (20) 

To avoid proliferation of indices, we will now write, e.g., as Di where 
i is now understood to run over just pixels, just antenna patterns, or both, 
according to context. The data is once again modelled as due to signal from 
the microwave background, Si and noise from the atmosphere and instrument 
iVj. The time stream, d ia is also split into signal s(8 ia ) and noise n ia . 

The two-point theoretical correlation function C£ = (SiSj) can be easily 
related to the intrinsic two-point correlation function C s (9 a i, 6^) = (s(9 a i)s(9bj))' 



(19) 
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C Z = lh) EE( 5 W s (^)>v* (2i) 

a b 



n 



By isotropy in the mean, C s (O a i, Obj) only depends on the angular distance, 

Qiajb, between 9 l a and 6{. Thus, the correlation function can be decomposed 
into Legendre polynomials and Eq. [ZT] can be rewritten: 

C^Y^C^, (22) 



where 



i 47r 



w m = (^) 2 J2J2 p d cos ( e ^ jb ))w ia w Jb - (23) 

a b 



n 



If the noise is white (uncorrelated from sample to sample) then the noise 
matrix 



and therefore 



NET% 

((«««]) = n ^^$ij (24) 



C£MW) = ^QE«#«- (25) 



A convenient normalization prescription for the weight vector is to set 
J2i w i/n to a constant because then the variance in the noise is the same 
for every antenna pattern. Taking that constant to be unity, we find the 
variance in the noise for each antenna pattern is simply NET 2 /t where t is 
the total observing time - the same formula as for "stare mode" . This weight 
vector normalization is equivalent to a window function normalization and 
is actually a very sensible one. Since the NET is the same for all antenna 
patterns, the window functions with the higher signal-to-noise ratios will 
have larger amplitudes. Thus plotting the window functions normalized in 
this way allows for quick graphical comparson of both the 1-space coverage 
and relative sensitivities of the various antenna patterns. 
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